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Problem  Statement: 

Matrix  eigenv'alue  [)robleins  play  a  significant  role  in  dynamic  analysis  of  structures.  The  natuial 
lrec[uencies  and  the  mode's  ol  a  Irc'e  vibrating  system  are  intiiTiatel\  related  to  the  eigenvalues 
and  the  eigenv’ectors  of  the  gein'ralized  system  ol  characteristic  equations  corresponding  to  that 
vibrating  sv’stem.  to  calculate  even  a.  few  ol  tlie  iiatiirai  frequencies  and  modes  of  a.  large  scale 
structure  is  computationallv’  intensive.  A  [recjuently  encountered  problem  in  structiiial  dynamics 
is  how  to  take  into  account  analysis  and  design  changes  introduced  after  the  structural  dynamic 
analysis  has  been  completed  and  tlie  natural  frequencies  and  modes  hav^e  been  computed.  The 
problem  is  how  the  natural  frequencies  and  modes  obtained  from  a  previous  analysis  can  be  used  to 
derive  the  response  of  the  new  modified  structure,  without  extensive  additional  computations.  In 
particular,  we  focus  on  the  structural  dynamic  re-analysis  problem  that  the  system  is  modified  in  a 
few  small  localized  regions  -  the  incremental  stiffness  matrix  and  the  incremental  mass  matrix  are 
low  rank  matrices,  but  t  heir  magnit  ude  need  not  l)e  small  relativ^e  to  the  original  matrices.  Dynamic 
solution  procedures  that  are  suitable  for  parallel  computer  environments  are  also  investigated. 

Summary  of  Result? 

In  this  project,  we  have  develoi)ed  and  implemented  a  number  of  parallel  algorithms  related  to 
matrix  factorization  and  generalized  eigenvalue  problems.  Furthermore,  a  Lanezos  based  method 
for  the  modified  eigenvalue  problem  has  beeji  developed.  The  following  summarizes  the  results  of 
this  project: 

1.  Development  and  imph'mentation  of  a  new  parallel  sparse  factorization  method(Ref.  1,3). 
Specific  contribution  includes  an  efficient  j)rocedure  for  the  solution  of  triangular  systems, 
which  is  generally  known  as  one  of  the  most  difficult  problems  for  the  parallelization  of  direct 
solution  methods.  The  parallel  procedure  for  the  solution  of  triangular  systems  is  particularly 
important  for  generalized  eigenvalue  problem  which  recjuires  solution  of  linear  systems  at  each 
Lanezos  iteration  step. 

2.  A  parallel  implementation  of  a  robust  Lanezos  algorithm  for  generalized  eigenvalue  problems 
(Ref.  1,10).  The  parallel  algorithm  has  been  implemented  in  a  finite  element  structural  anal- 
V'sis  program.  To  ensure  the  robustness  and  accuracy,  the  parallel  procedure  implemented 
includes  spectral  sliilts,  partial  and  selective  re-orthogonalizations.  We  have  also  tested  vaii- 
ous  approaches  and  developed  a  rol)ust  implementation  so  that  tlie  finite  element  matrices  can 
be  in  the  form  of  ('irlu'r  eleiiu'iit  matrices  or  assembled  system  matrices.  Tradeoffs  between 
spectral  shifts,  factorization  and  re-ortliogonalization  have  been  studied. 
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i.  The  developnient  of  a  new  niotliod  lor  scdviiig  perturl)ed  syninietric  eigenvalue  problems 
(Ref.  3.5.9).  This  met  hod  is  based  on  an  extension  of  the  block  Lanezos  method  for  standard 
and  generalized  eigcMivahie  j)roblems.  This  new  method  is  capable  of  updating  the  modified 
eigenvalues  and  eigenvectors  even  when  the  modification  made  to  the  structures  may  cause 
significant  changes  to  the  mode  shapes  of  the  system. 

In  addition,  we  ha\'e  also  studied  the  following  ])roblems  related  to  eigenvalue  problems  and  the 
Lanezos  and  conjugate  gradient  methods  for  structural  analysis: 

1.  A  study  of  the  solution  methods  for  the  minimization  of  the  generalized  Rayleigh  Quotient 
subject  to  linear  and  quadratic  constraints  (Ref.  7).  Applications  of  the  methods  for  the 
evaluation  of  thin  sliell  finite  elements  are  illustrated. 

2.  Development  of  a  robust  and  sta!)le  incomplete  Cholesky  (IC)  factorization  (Ref.  2,8).  In 
particular,  we  have  examined  in  detail  the  tradeoffs  between  the  storage  requirement  and  the 
effectiveness  (in  terms  of  the  total  CPU  time  for  the  solution)  of  the  IC  preconditioners  for 
the  conjugate  gradient  metliod.  A  technique  lias  been  developed  to  manage  both  the  storage 
space  and  to  optimize  t  he  “goodness''  of  tlu^  preconditioning  matrix. 

3.  Development  of  a  simi)le  but  effective  procedure  to  enhance  the  block  conjugate  gradient 
method  for  the  solution  of  a  linear  system  of  equations  with  multiple  right  hand  sides  (Ref. 
2,6).  This  development  is  important  since  engineering  structures  are  often  designed  to  resist 
multiple  load  cases. 

4.  A  studv  of  the  use  of  the  conjugate  gradient  method  for  nonlinear  analysis  of  structures 
(including  snap  thi’ough  and  simple  l)ifurcation  of  structures)  (Ref.  2).  Based  on  the  rela¬ 
tionship  between  the  conjugate  gradient  and  t  he  Lanezos  methods,  we  have  examined  how  to 
take  acl vantage  of  the  information  (such  as  the  eigenvalues,  the  residual  and  search  direction 
vectors)  in  soh’ing  nonlinear  j>rol)lems.  For  example.  l)y  simply  keeping  track  of  the  inertia 
of  the  (tangential)  stiffness  matrix,  we  can  determine  easily  the  secondary  load-displacement 
path  related  to  simple  nonlinear  liifurcation  of  structures. 

5.  A  comparison  study  on  a  vai*iety  of  preconditioners  (including  polynomial  preconditioners, 
incomplete  Cholesky  factors,  domain  decomposition  and  a  combination  of  the  preconditioning 
methods)  for  conjugate  gradient  methods  on  a  distributed  memory  parallel  computer  (Ref. 
2).  Both  linear  and  nonlinear  structural  problems  have  been  examined  in  this  study. 

All  parallel  solution  procedures  have  been  implemented  and  tested  on  an  Intel  iS60  parallel  com¬ 
puter. 
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